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METHOD AND DEVICE GENERATING DATA INPUT TO AN 
ELECTROMAGNETIC FIELD INTENSITY CALCULATING DEVICE 

Background of the Invention 
5 Field of the Invention 

The present invention relates to a method 
calculating the intensity of an electromagnetic field 
emitted from an electric circuit device, and more 
particularly, to a method and device generating data 
10 input to an electromagnetic field intensity calculating 
device which calculates the intensity of an 
electromagnetic field emitted from an electric circuit 
device by using a moment method. 

15 Description of the Related Art 

The electromagnetic field intensity of an 
electromagnetic wave emitted from an object having an 
arbitrary shape, for example, an electric circuit device 
can be easily calculated if an electric current which 

2 0 flows in each part of the object is known. Theoretically, 
the value of an electric current is given by solving 
the Maxwell's electromagnetic field equations under a 
given boundary condition. However, there is no 
theoretical solution to an object having an arbitrary 

25 shape under a complex condition. 
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Accordingly, all of electric current calculation 
methods used in a current electromagnetic field 
intensity calculating device are approximation methods. 
Three methods such as an infinitesimal loop antenna 
5 approximation method, a distributed constant line 
approximation method, and a moment method are known as 
H the approximation solutions. 

Q The infinitesimal loop antenna approximation 

%. method is a method with which a wire interconnecting 

^1 10 a wave source circuit and a load circuit is handled as 

a loop antenna, and an electric current on the loop is 
obtained with a calculation of a lumped constant circuit 
by assuming that the electric current on the loop is 
flat. Although the calculation of this method is the 

15 simplest, its accuracy is low under a condition that 
the size of a loop cannot be ignored in comparison with 
the wavelength of an electromagnetic wave. 

The distributed constant line approximation 
method is a method with which an electric current is 

2 0 obtained by applying the distributed constant line 
equations to an object that can be approximated to a 
one-dimensional structure. Its calculation is 
relatively easy, and an electric current can be quickly 
obtained with high accuracy. However, for an object to 

25 be analyzed, which cannot be approximated to a 
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one-dimensional structure, by way of example, an object 
having an arbitrary shape such as a printed circuit board, 
a metal cabinet, etc., an electric current cannot be 
obtained. 

5 In the meantime, the moment method is one solution 

of integral equations derived from the Maxwell's 
electromagnetic field equations, and can handle a 
three-dimensional object having an arbitrary shape. 
With this method, an object is partitioned into small 

10 elements (patches or meshes), and an electric current 
is calculated. 

If the intensity of an electromagnetic field 
emitted from an electric circuit device is calculated 
with the moment method, it can be calculated with high 

15 accuracy by combining all of constituent elements such 
as a printed circuit board, a cable, a lead, a metal 
cabinet, etc. within the electric circuit device, and 
by accurately putting the device into a model. 

Three-dimensional data and data of wires are 

20 captured from a CAD as data of a printed circuit board, 
a cable, a lead, etc., so that an analysis model can 
be generated with relative ease. However, for 
three-dimensional structure data such as the metal 
cabinet of an electric circuit device, etc., ameshmodel 

25 must be generated by partitioning the data into meshes 



(patches) which are suitable for the moment method after 
capturing the data from a CAD. At this time, the model 
must be generated accurately and efficiently. 

Conventionally, a metal cabinet of an electric 
circuit device, etc. is conventionally partitioned into 
meshes (patches) shaped like triangles each of which 
is configured by 3 points, as an automatic partitioning 
method, even if the metal cabinet, etc. is partitioned 
into meshes which are suitable for the moment method. 
As will be described later, a mesh partitioning method 
that is most suitable for the moment method is a method 
partitioning into quadrilaterals. Conventionally, 
however, it is difficult to automatically partition a 
metal cabinet, etc. into quadrilateral meshes each of 
which is configured by four points. 

That is, conventionally, for a portion that cannot 
be partitioned as a quadrilateral, or a portion having 
a hole, CAD data must be manually changed so that a CAD 
model must be modified to a form suitable for 
partitioning into quadrilateral meshes. This requires 
much time and labor. 

Figs. 1, 2A, and 2B explain the influences that 
different mesh shapes exert on an electric current 
calculation using the moment method. Fig. 1 explains 
the difference between electric current calculations 



within meshes (patches) . In this figure, electric 
currents flowing from vertexes to the directions 
(indicated by arrows) of opposite sides are calculated 
for a triangle mesh, whereas electric currents flowing 
to the directions of opposite sides are calculated for 
a quadrilateral mesh. 

Figs . 2A and 2B explain the directions of electric 
currents as a whole in the case where such meshes are 
combined. If partitioning is made into triangle meshes 
as shown in Fig. 2A, the flow of an electric current 
becomes zigzag, namely, uneven, and the length of a path 
on which the electric current flows becomes longer. As 
a result, a propagation delay of the electric current 
occurs, and the accuracy of an analysis made with the 
moment method is degraded. In the meantime, if 
partitioning is made into quadrilateral meshes, an 
electric current smoothly flows, and the accuracy of 
an analysis made with the moment method is improved. 

Siommary of the Invention 

An object of the present invention is to improve 
the analysis accuracy of an electromagnetic field 
intensity calculation as a whole, and to increase an 
efficiency of the process of the electromagnetic field 
intensity calculation in view of the above described 
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problem by automatically partitioning a metal cabinet 
of an electric circuit device, etc. into quadrilateral 
meshes which are suitable for an electric current 
calculation using the moment method, and by providing 
5 the partitioned data to an electromagnetic field 
intensity calculating device. 

A method generating data input to an 
electromagnetic field intensity calculating device 
according to the present invention, which is a method 

10 generating data input to an electromagnetic field 
intensity calculating device that calculates the 
intensity of an electromagnetic field emitted from an 
electric circuit device having a metal cabinet, extracts 
surface data of the metal cabinet from three-dimensional 

15 data of the electric circuit device, partitions the 
surface corresponding to the surface data into 
quadrilateral meshes, and outputs the data partitioned 
into meshes to the electromagnetic field intensity 
calculating device. 

20 According to the present invention, a metal 

cabinet model which is suitable for a calculation of 
the moment method, namely, a model partitioned into 
quadrilateral meshes can be automatically and quickly 
generated from the three-dimensional structure data of 

25 a CAD system, and such a model is output to an 
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electromagnetic field intensity calculating device, 
thereby accurately calculating the intensity of an 
electromagnetic field emitted from the whole of an 
electric circuit device including its metal cabinet. 
5 This can greatly contribute to an improvement in the 
practicality of the electromagnetic field intensity 
calculating device. 

Brief Description of the Drawings 

10 Fig. 1 explains the difference between electric 

current calculations depending on mesh shapes; 

Figs . 2A and 2B explain the difference between the 
directions of electric currents flowing in triangle and 
quadrilateral meshes; 
15 Fig. 3 is a block diagram showing the principle 

of the present invention; 

Fig. 4 shows the relationship between an input 
data generating device according to the present 
invention and an electromagnetic field intensity 
20 calculating device; 

Fig. 5 is a block diagram showing the 
configuration of the input data generating device; 

Fig. 6 explains the method representing a curved 
surface in a CAD data file; 
25 Fig. 7 shows the data structure of a CAD data file; 
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Fig. 8 explains the partitioning of a metal 
cabinet into quadrilateral meshes on the metal cabinet 
surface; 

Fig. 9 explains a same surface existing on the 

front or the back of a metal cabinets- 
Fig . 10 explains the method extracting four 

vertexes to approximate each surface to a quadrilateral 

from control points representing each surface 

configuring the same surface; 

Fig. 11 explains the method equally partitioning 

a flat or a curved surface, which has the four points 

obtained in Fig. 10 as vertexes, for opposite sides; 

Fig. 12 explains the method sharing partitioning 

points on a connecting side of contiguous surfaces; 

Fig. 13 explains the method repartitioning meshes 

in the case where surfaces of different materials are 

superposed; 

Fig. 14 is a flowchart showing the details of a 

same surface extraction process; 

Fig. 15 explains a same surface as a curved 
surface; 

Fig. 16 explains the progression (No. 1) of the 
process for extracting a same surface from surface data 
of a metal plate having a curved surface; 

Fig. 17 explains the progression (No. 2) of the 



process for extracting a same surface from surface data 
of a metal plate having a curved surface; 

Fig. 18 explains the progression (No. 3) of the 
process for extracting a same surface from surface data 
of a metal plate having a curved surface; 

Fig. 19 explains the progression (No. 4) of the 
process for extracting a same surface from surface data 
of a metal plate having a curved surface; 

Fig. 20 explains a same surface of a metal plate 
having a flat surface; 

Fig. 21 explains the progression (No. 1) of a same 
surface extraction process for a metal plate having a 
flat surface (No. 1); 

Fig. 22 explains the progression (No. 2) of the 
same surface extraction process for a metal plate having 
a flat surface; 

Fig. 23 explains the progression (No. 3) of the 
same surface extraction process for a metal plate having 
a flat surface; 

Fig. 24 explains the progression (No. 4) of the 
same surface extraction process for a metal plate having 
a flat surface; 

Fig. 25 is a flowchart showing the details of a 
4-point block surface extraction/mesh partitioning 
process; 
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Fig. 26 exemplifies the extraction of a 4-point 
block surface (No. 1); 

Fig. 27 exemplifies the extraction of the 4-point 
block surface (No. 2) ; 
5 Fig. 28 exemplifies the extraction of the 4-point 

block surface (No. 3) ; 

Fig. 29 exemplifies the extraction of the 4-point 
^ block surface (No. 4} ; 

Fig. 30 explains the method representing a curved 
^ 10 surface in data partitioned into meshes; 

„ Fig. 31 shows the structure of data partitioned 

'7 into meshes; 

I'jf Fig. 32 is a flowchart showing the details of the 

- process for a point shared by contiguous surfaces; 

15 Fig. 33 shows a specific example of the process 

for a point shared by contiguous surfaces; 

Fig. 34 is a flowchart showing the details of the 
process for aligning a mesh shape; 

Fig. 35 exemplifies mesh shape alignment; and 
20 Fig. 3 6 explains the loading of a program 

according to a preferred embodiment into a computer. 



Description of the Preferred Embodiments 

Fig. 3 is a block diagram showing the principle 
25 of the present invention. This figure is a block diagram 
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showing the functions of a method generating data input 
to an electromagnetic field intensity calculating 
device which calculates the intensity of an 
electromagnetic field emitted from an electric circuit 
5 device having a metal cabinet. 

In this figure, with the input data generating 
method according to the present invention, surface data 
of a metal cabinet is extracted from three-dimensional 
data of an electric circuit device in 1, the surface 

10 corresponding to the surface data is partitioned into 
quadrilateral meshes in 2, and the data partitioned into 
meshes is output to an electromagnetic field intensity 
calculating device in 3. 

In a preferred embodiment of the present invention, 

15 if surface data of a metal cabinet is composed of a 
plurality of pieces of surface data including the front 
and the back of a metal plate which configures the metal 
cabinet, data of a same surface composed of a plurality 
of surfaces that exist on the front or the back of the 

20 metal plate is extracted from the plurality of pieces 
of surface data, and each of the surfaces configuring 
the same surface can be also partitioned into 
quadrilateral meshes. 

Additionally, control points corresponding to the 

25 vertexes of a quadrilateral can be extracted by using 
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the data of control points as the surface data of each 
of surfaces which configure the same surface, when each 
of the surfaces is approximated to the quadrilateral. 
Furthermore, the control points corresponding to the 
vertexes of the quadrilateral which approximates each 
of the surfaces are recognized as 4 vertexes, and a flat 
or a curved surface determined by the data of the control 
points is equally partitioned for respective pairs of 
opposite sides so that the surface can be partitioned 
into quadrilateral meshes. 

Furthermore, in the preferred embodiment 
according to the present invention, if surface data of 
a metal cabinet is data corresponding to a plurality 
of surfaces, each of the plurality of surfaces is 
partitioned into quadrilateral meshes, and a comparison 
is made between the coordinates of partitioning points 
on two sides that can possibly be a side shared by 
contiguous surfaces among the plurality of surfaces. 
If the coordinates of the partitioning points are 
determined to match within a preset tolerance, the data 
of the partitioning points are recognized to be data 
shared by the contiguous surfaces. The data which is 
partitioned into meshes and includes the shared data 
can be output to the electromagnetic field intensity 
calculating device. 



Furthermore, in the preferred embodiment, if a 
surface whose material is different from a metal cabinet, 
and which is superposed on a surface corresponding to 
surface data of the metal cabinet and has an area smaller 
than the surface corresponding to the surface data, 
quadrilateral meshes are repartitioned by being aligned 
with the shape of the surface having the smaller area, 
after the surface corresponding to the surface data is 
partitioned into the quadrilateral meshes. The 
repartitioned data can be output to the electromagnetic 
field intensity calculating device. 

Used as a storage medium according to the present 
invention is a computer-readable storage medium used 
by a computer generating data input to an 
electromagnetic field intensity calculating device 
that calculates the intensity of an electromagnetic 
field emitted from an electric circuit device having 
a metal cabinet, on which is recorded a program for 
causing the computer to execute a process, the process 
comprising: receiving specification of one or more 
surfaces among surfaces configuring the metal cabinet 
from a user; extracting surface data of the specified 
surface from the three-dimensional data of the electric 
circuit device; partitioning the surface corresponding 
to the surface data into quadrilateral meshes; and 
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outputting the data partitioned into meshes to the 
electromagnetic field intensity calculating device. 

A device according to the present invention, which 
generates data input to an electromagnetic field 
5 intensity calculating device, comprises a surface data 
extracting unit extracting surface data of a metal 
cabinet from three-dimensional data of an electric 
circuit device, a mesh partitioning unit partitioning 
the surface corresponding to the surface data into 

10 quadrilateral meshes, and a generated data outputting 
unit outputting the data partitioned into meshes to the 
electromagnetic field intensity calculating device. 

As described above, according to the present 
invention, a surface corresponding to the surface data 

15 of a metal cabinet is partitioned into quadrilateral 
meshes which are suitable for a calculation of the moment 
method in order to generate data input to an 
electromagnetic field intensity calculating device 
which calculates the intensity of an electromagnetic 

20 field emitted from an electric circuit device having 
the metal cabinet . 

Fig. 4 explains the relationship between an 
electromagnetic field intensity calculating device and 
an input data generating device according to the present 

25 invention, which generates data input to the 



electromagnetic field intensity calculating device. 
The output of the input data generating device 11 is 
provided to the electromagnetic field intensity 
calculating device 12, which then calculates the 
intensity of an electromagnetic field emitted, for 
example, from an electric circuit device. 

A preferred embodiment of the present invention 
is explained by assuming that the input data generating 
device 11 partitions, for example, a metal cabinet of 
an electric circuit device into quadrilateral meshes 
which are suitable for a calculation of the moment method, 
and provides the data partitioned into meshes to the 
electromagnetic field intensity calculating device 12. 
However, the input data generating device 11 may 
generate input data corresponding not only to the metal 
cabinet, but also to constituent elements such as a 
printed circuit board, a cable, a lead, etc. within the 
electric circuit device. 

The electromagnetic field intensity calculating 
device 12 is configured by an impedance calculating unit 
13 calculating the self-impedance of each partitioned 
mesh (patch) or the mutual impedance between meshes, 
etc. in correspondence with an output of data 
partitioned into meshes, for example, for a metal 
cabinet; a current distribution calculating unit 14 
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calculating an electric current in each portion of an 
electric circuit device by using a result of the 
calculation made by the impedance calculating unit 13, 
and the voltage of a wave source, which is provided, 
5 for example, from a user; and an electromagnetic field 
calculating unit 15 calculating with an obtained current 
distribution an electric or a magnetic field at a point 
apart, for example, from the electric circuit device 
by a predetermined distance. 
10 Fig. 5 is a block diagram showing the 

configuration of a device generating data input to an 
electromagnetic field intensity calculating device, 
according to a preferred embodiment of the present 
invention. 

15 In this figure, three-dimensional structure data 

of a metal cabinet is read by a three-dimensional 
structure data reading unit 21 from a CAD data file 20, 
data of all surfaces are sequentially stored in a memory 
unit 22 as surface data, and the shape of each of the 

20 surfaces is displayed by a display unit 23 . For example, 
a user selects a surface to be partitioned into meshes 
with a click via a keyboard/mouse input unit 24. The 
selected surface is then provided to a same surface 
extracting unit 25. 

25 The same surface extracting unit 25 extracts a 
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same surface if data of each surface stored in the CAD 
data file 20 is block data. The block data includes the 
data of the front, the back and the sides of each metal 
plate that configures a metal cabinet. In this case, 
5 it is sufficient to partition only either of the front 
and the back of the metal plate into meshes, and it is 
unnecessary to partition the opposite surface in the 
direction of the thickness of the metal plate, and the 
sides. Therefore, a plurality of surfaces configuring 

10 the front and the back are extracted as a same surface. 
If data stored in the CAD data file 20 is surface data 
that ignores the weight of a metal plate, this extraction 
process becomes unnecessary. The same surface 
extraction process will be described in detail later. 

15 A result of the process performed by the same 

surface extracting unit 25 is provided to a 4-point block 
surface extracting/mesh partitioning unit 26. 
According to this preferred embodiment, one surface is 
approximated to a quadrilateral as much as possible, 

20 and the quadrilateral is equally partitioned into 
quadrilateral meshes for opposite sides. With the 
4-point block surface extraction, 4 vertexes are 
extracted so that the shape of the surface before being 
partitioned into meshes becomes as close to a 

25 quadrilateral as possible. An extracted 4-point block 



18 



surface is equally partitioned for opposite sides by 
using a partitioning length input via the keyboard/mouse 
inputting unit 24, for example, from a user. 

Then, a contiguous surface point sharing unit 27 
5 makes a comparison between the coordinates of 
partitioning points corresponding to a side shared by 
contiguous surfaces. Partitioning points that are 
Ci determined to match, for example, within a preset 

%j tolerance are recognized as data shared by the 

Si; 

10 contiguous surfaces, and meshes are generated for the 

contiguous surfaces connected. 
!p' Next, a mesh shape aligning unit 28 repartitions 

the meshes in correspondence with a superposed portion 
O of a surface of a different material. Meshes are 

Hi 

15 repartitioned to align the shape of a hole, for example, 
if a metal plate has the hole. 

Generated mesh data is stored in a mesh data file 
30 by a mesh data writing unit 29. Contents of this file 
are output as data input to the electromagnetic field 

20 intensity calculating device, and used to analyze a 
radio wave leakage from a chink of a cabinet, or a shield. 
The process from the three-dimensional structure data 
reading unit 21 to the mesh data writing unit 29 is 
performed by using the working area of the memory unit 

25 22 on demand. 



The operations of the input data generating device 
shown in Fig. 5 are further described with reference 
to Figs. 6 to 13. Figs. 6 and 7 exemplify data stored 
in the CAD data file 20. Fig. 6 exemplifies a surface 
configuring a metal cabinet, normally, a curved surf ace . 
The curved surface is represented as a B-spline curved 
surface by control points of the curved surface. Each 
of the control points is represented by a parameter u 
in a T direction, a parameter v in an S direction, and 
a weight. Namely, a reference point having x, y, and 
z coordinates exists, and each of the control points 
has (displacement) coordinates u and v on each surface 
from the reference point. 

Fig. 7 exemplifies the structure of CAD data 
corresponding to the data shown in Fig. 6. This data 
includes the surface number of the B-spline curved 
surface, spatial coordinate values X, Y, and Z of each 
of the control points, and the like. 

Fig. 8 explains the mesh partitioning performed 
by the 4-point block surface extracting/mesh 
partitioning unit 26. Surface data of a metal cabinet 
is extracted from the three-dimensional structure data 
stored in the CAD data file 20, and the surface 
corresponding to the surface data is partitioned into 
quadrilateral meshes on the cabinet surface. A square 
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or a rectangle is most suitable as the quadrilateral. 
However, a trapezoid, a rhombus, etc. may be used. 

Fig. 9 explains the process performed by the same 
surface extracting unit 25. A user specifies (selects) 
5 any of a plurality of surfaces configuring a metal 
cabinet via the keyboard/mouse input unit 24 shown in 
Fig. 5. The same surface extracting unit 25 extracts 
the data of the same surface existing on the front or 
the back of the metal cabinet depending on whether the 
10 specified surface exists either on the front or on the 
back. 

Figs. 10 and 11 explain the process performed by 
the 4-point block surface extracting/mesh partitioning 
unit 26. In Fig. 10, for example, control points that 
15 become four vertexes for approximating each surface to 
a quadrilateral are extracted from the control points 
that represent each of the surfaces which are extracted 
in Fig. 9 and configure the same surface. 

In Fig. 11, the 4 control points extracted as shown 
20 in Fig. 10 are recognized as 4 vertexes of the 
quadrilateral, and the process for equally partitioning 
the flat or the curved surface configured by these 4 
vertexes for opposite sides is performed. 

Fig. 12 explains the process performed by the 
25 contiguous surface point sharing unit 27. If coordinate 



values of mesh partitioning points on sides that are 
considered to be a side shared by contiguous surfaces 
to be connected are determined to be the same within 
a tolerance, these surfaces are recognized to be 
connected by this side, and the partitioning points are 
shared- 
Fig. 13 explains the process performed by the mesh 
shape aligning unit 28 when two surfaces of different 
materials are superposed. In this figure, meshes are 
repartitioned by being aligned with the shape of a 
surface having a smaller area when the surfaces of 
different materials are superposed. For example, the 
upper left square mesh among 9 large square meshes is 
to be repartitioned into 3 rectangular meshes on the 
right side of Fig. 13. 

Fig. 14 is a flowchart showing the details of the 
same surface extraction process performed by the same 
surface extracting unit 25. Once the process starts in 
this figure, a surface configuring a metal cabinet is 
clicked with a mouse, for example, by a user as a 
reference surface in step SI. In step S2, the reference 
surface is selected. Then, in step S3, surface data is 
output from the memory unit 22 shown in Fig. 5. In step 
S4, the reference surface, namely, the initial process 
target surface, or a surface contiguous to the process 
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target surface is extracted. 

Here, the reference surface is first extracted. 

In step S5, the barycenter of the reference surface, 

and a normal in the outer direction of the metal plate 

5 are calculated. Then, in step S6, a vector in the 

direction opposite to the normal is calculated. In step 

S7, a surface which intersects in the opposite direction 

^. is searched by using the surface data. 

Such a process is explained with reference to Figs . 

IS! 10 15 to 24. Fig. 15 exemplifies a curved surface of a metal 

cabinet, whereas Fig. 20 exemplifies a flat surface of 

a metal cabinet. When a user clicks an upper surface 

t'y on the left side with a mouse as a reference surface 

~~ in the example shown in Fig. 15, a normal vector for 

PJ 

15 the selected reference surface A, and a vector in the 
opposite direction are calculated as shown in Fig. 16. 

A surface intersecting the vector in the opposite 
direction is searched in step S7 of Fig. 14. In Fig. 
16, the back of the reference surface A is searched. 

20 In step S8 of Fig. 14, it is determined whether or not 
a distance to the intersecting surface is within a preset 
tolerance. Here, the distance is the thickness of the 
metal plate. If it is determined that the distance is 
within the tolerance, it is further determined in step 

25 S9 whether or not the intersecting surface is a surface 
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which has been determined as the same surface. 

In Fig. 15, the surfaces to be determined as the 
same surface is the reference surface, and the surface 
right to the reference surface, namely, the upper 
5 surface on the right side. Since the process has just 
started at this time, the intersecting surface is 
determined not as the same surface. Then, in step SIO, 
it is determined whether or not the intersecting surface 
is a surface which has been removed. Here, the surface 

10 is also determined not as the same surface. In step Sll, 
a removal flag is set for this surface, and the back 
of the reference surface A is removed. 

In step S12, a same surface flag is set for the 
process target surface, here, the reference surface A. 

15 In step S13, this surface, here, the reference surface 
A is highlighted on the display. In step S14, it is 
determined whether or not the next surface contiguous 
to the reference surface A exists. 

A surface B contiguous to the reference surface 

20 A exists as shown in Fig. 17, and the operations in and 
after step S4 are performed for this contiguous surface 
recognized as a process target surface. A surface which 
intersects the vector in the direction opposite to the 
normal vector of the process target surface B is searched 

25 in a similar manner. After the operations in steps S8 



to SIO are performed, the removal flag is set also for 
the intersecting surface in step Sll, and the same 
surface flag is set for the process target surface, 
namely, the surface B in step S12. In step S13, the 
process target surface is highlighted. 

Then, in step S14, it is determined whether or not 
the next contiguous surface exists. Here, a contiguous 
surface is not necessarily contiguous to an immediately 
preceding process target surface, and may be considered 
as a surface yet to be processed. For example, a side 
C is selected as shown in Fig. 18, and the operations 
in and after step S4 are performed. 

In a similar manner as in the above described case, 
a surface intersecting an opposite direction vector is 
searched in step S7 . Here, the surface A as the reference 
surface selected by the user is found as the intersecting 
surface. Since the same surface flag is already set for 
the surface A, this surface is determined to be a surface 
that has determined as the same surface in step S9. In 
step S15, the removal flag is set for the process target 
surface, namely, the side C itself. The side C is then 
removed. 

In step S14, a side D is selected as a surface yet 
to be processed as shown in Fig. 19, and the operations 
in and after step S4 are performed. A surface that is 
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searched as a surface intersecting an opposite direction 
vector in step S7 is a surface that has been removed 
as explained with reference to Fig. 17. In step SIO, 
it is determined that the searched surface is a surface 
5 which has been removed. Then, in step S15, the process 
target surface, namely, the side D is also removed. Here, 
the process is terminated for all of the surfaces. 
.Z.' Surfaces determined to be the same surface as the 

"J' reference surface are two upper surfaces as shown on 

\\ 

If: 10 the right side of Fig. 15. 

J| Figs. 21 to 24 explain the process for a flat 

surface shown in Fig. 20. A difference from the process 

^' shown in Figs. 16 to 19 is the process shown in Figs. 

Ill 

^ 23 and 24. In Fig. 23, the process is performed for a 

O 

fy 15 side C on the left side, which is recognized as a process 

target surface. In step S8 of Fig. 14, the distance to 
the surface intersecting an opposite direction vector, 
here, a surface D is determined to be longer than a 
tolerance. Therefore, the process target surface, 

20 namely, the surface C itself is removed in step S15. 

If the operation of step S8 in Fig. 14 is not 
performed, or if the distance to the intersecting 
surface is not longer than the thickness of the metal 
plate and not shorter than the width (the distance 

25 between the sides) , the same surface flag is set for 
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the process target surface, here, the side C in step 
S12 after the operations in steps S9 to Sll are performed. 
Accordingly, the side C is determined to be the same 
as the reference surface A. 
5 In Fig. 24, the process is performed for the side 

D on the right side, which is recognized as a process 
target surface. Similar to the side C, the distance to 
the intersecting surface is determined to be longer than 
the tolerance in step S8. As a result, the removal flag 

10 is set for the process target surface, namely, the side 
D in step S15. 

If a metal cabinet plate has a curved surface as 
shown in Fig. 15, the same surface extraction process 
explained in this preferred embodiment does not always 

15 yield a proper result. For example, if not the surface 
B but the side C is determined as a process target surface, 
and the vector in the direction opposite to the surface 
C intersects not the reference surface Abut the surface 
B after the same surface flag is set for the reference 

20 surf ace A in Fig. 16, also the side surf ace C can possibly 
be determined as the same surface. Therefore, it should 
be noted that a proper result cannot be obtained 
depending on the shape of a curved surface, the 
arrangement of each surface configuring the same surface, 

25 the setting of a tolerance for the distance to an 
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intersecting surface in step S8 of Fig. 14, or the like. 

Fig. 25 is a flowchart showing the details of the 
process performed by the 4-point block surface 
extracting/mesh partitioning unit 2 6 shown in Fig. 5. 
5 Once the process starts in this figure, one of a 
plurality of surfaces which configure a same surface 
is extracted from the same surface in step S22 by using 
a mesh partitioning length input from a user in step 
S20, and same surface data output in step S21. In step 

10 S23, it is determined whether or not the process is 
completed for all of the plurality of surfaces that 
configure the same surface. If the process is not 
completed, distances from the barycenter of the 
extracted surface to a plurality of control points 

15 representing the surface are calculated for the 
extracted surface in step S24. Then, in step S25, a 
control point whose distance is the longest is 
determined to be a reference point - 

Subsequent steps S2 6 to S30 correspond to the 

20 4-point block surface extraction process, which is 
described below by citing the examples shown in Figs. 
26 to 29. 

In step S2 6, all of control points within a range, 
which centers the barycenter, of 180 degrees from a 
25 reference point, a point A in Fig. 26 are sequentially 
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extracted. Then, in step S27, the area of a triangle 
which is configured by each of the extracted control 
points, the barycenter, and the reference point is 
calculated in step S27. In step S28, it is determined 

5 whether or not all of the control points within the range 
of 180 degrees have been extracted. If all of the control 
points have not been extracted, the operations in and 
after step S26 are repeated. 

If it is determined that all of the control points 

10 within the range of 180 degrees have been extracted, 
the control point corresponding to the triangle having 
the largest area is recognized as the next reference 
point in step S29. Then, in step S30, it is determined 
whether or not 4 points of a 4-point block surface, 

15 namely, the 4 vertexes of the quadrilateral have been 
extracted as reference points. If these points have not 
been extracted, the operations in and after step S26 
are repeated. 

In Fig. 26, points B, C, and D are sequentially 

20 extracted with the loop in steps S26 to S28 as all of 
the control points within the range of 180 degrees from 
the point which is recognized as the reference point 
in step S25, namely, the point A whose distance from 
the barycenter is determined to be the longest. In step 

25 S27, the areas of respective triangles are calculated. 



In step S29, the control point corresponding to the 
triangle having the largest area, here, the point B is 
recognized as the next reference point. Assuming that 
the barycenter is a point 0 in Fig. 26, the area of a 
triangle AOB is supposed to be larger than that of a 
triangle AOC. 

Fig. 27 explains the operations in steps S26 to 
S30 by recognizing the point B as a reference point. 
Points C, D, and E are extracted as control points within 
the range, which centers the barycenter O, of 180 degrees 
from the reference point B, and the control point E 
corresponding to the triangle having the largest area 
is recognized as the next reference point. 

In Fig. 28, control points within the range, which 
centers the barycenter, of 180 degrees from the 
reference point E are extracted. If control points 
within the range, which centers the barycenter, of 180 
degrees from a reference point as described above, a 
control point that has been recognized as a reference 
point is not extracted. Accordingly, a point F is 
recognized as the next reference point in step S29. In 
step S30, it is determined that 4 points have already 
been extracted, and a partitioning number of 4 sides 
is calculated in step S31 . 

The partitioning number is calculated in 
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correspondence with the mesh partitioning length input 
from a user in step S20, namely, the length of the side 
of a mesh. However, if a remainder occurs when the length 
of each side of a quadrilateral as a 4-point block is 
5 divided by a specified partitioning length, 1 is added 
to the partitioning number. 

Note that the mesh partitioning length input from 
a user is set to, for example, on the order of 1/20 of 
the wavelength of a maximum frequency in order to suit 

10 a calculation of the moment method. Additionally, if 
the partitioning numbers of opposite sides differ, also 
a shorter side is partitioned by the partitioning niimber 
of a longer side. 

Then, in step S32, the equal partitioning for 

15 opposite sides, which is explained with reference to 
Fig. 11, is performed. In Fig. 11, one pair of opposite 
sides is partitioned into 5 meshes, whereas the other 
is partitioned into 3 meshes. Then, in step S33, the 
mesh data are stored. After the mesh data are stored, 

20 one remaining surface is extracted from the same surface 
data in step S22. The operations in and after step S23 
are then performed. When a remaining surface cannot be 
extracted in step S22, namely, when it is determined 
that the process is completed for all of the surfaces 

25 in step S23, the process is terminated. 



Figs. 30 and 31 explain the structure of mesh data 
stored in the mesh data file 30 shown in Fig. 5. Fig. 
30 exemplifies quadrilateral meshes (patches) 
corresponding to the curved surface of a metal cabinet, 
5 and each of the meshes is represented by points 
corresponding to mesh partitioning points. 

Fig. 31 exemplifies the structure of mesh data. 
The mesh data is composed of coordinate specification 
data of vertexes of a polygon, and specification data 

10 of polygon configuring points. The coordinate 
specification data is composed of a point number for 
each of the points in Fig. 30, and x, y, and z coordinate 
values. The specification data of polygon configuring 
points is composed of a patch number for one polygon 

15 (quadrilateral mesh/patch) , and point numbers of points 
which are the 4 vertexes of the patch. 

Fig. 32 is a flowchart showing the details of the 
process for a mesh partitioning point shared by 
contiguous surfaces, namely, the process performed by 

20 the contiguous surface point sharing unit 27 shown in 
Fig. 5. Once the process starts in this figure, 
coordinate values of a mesh partitioning point are 
extracted from mesh data in step S41. In step S42, it 
is determined whether or not the process is completed 

25 for the coordinate values of all of partitioning points . 
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If the process is not completed, an interval between 
the coordinate values of the extracted partitioning 
point and a partitioning point which configures a 
contiguous surface is calculated in step S43. In step 
S44, it is determined whether or not the interval is 
within a tolerance. If the interval is larger than the 
tolerance, the operations in and after the calculation 
of the interval between the coordinate values of an 
extracted partitioning point and another partitioning 
point of the contiguous surface are repeated. In the 
calculation process of the interval between coordinate 
values, the x, y, and z coordinate values of a point 
corresponding to each partitioning point are already 
obtained as explained with reference to Fig. 31, and 
the interval is calculated by using these coordinate 
values . 

If it is determined that the interval between the 
coordinate values is within the tolerance in step S44, 
the two partitioning points are set as a shared point 
in step S45, and stored as mesh data. The process then 
goes back to step S41 where the coordinate values of 
another mesh partitioning point are extracted, and the 
operations in and after step S42 are repeated. If there 
are no coordinate values of a mesh partitioning point 
to be extracted, it is determined that the process is 
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completed for all of the partitioning points in step 
S42. Here, the point sharing process is terminated. 

Fig. 33 exemplifies a result of the point sharing 
process. If it is determined that the intervals between 
5 partitioning points of two sides are within a tolerance 
between the sides of two surfaces as shown in Fig. 33, 
these two sides are determined to be connecting sides 
of the two surfaces. The partitioning points of the two 
sides are therefore set as shared points. 
'\ 10 The process for sequentially extracting 4 

I vertexes of a 4-point block surface as reference points 

is earlier explained with reference to Fig. 25. If a 
shared point exists between contiguous 4-point block 
surfaces, this point is preferentially extracted as the 
IV 15 next reference point. 

Fig. 34 is a flowchart showing the details of the 
process performed by the mesh shape aligning unit 28 
shown in Fig. 5. With this process, if surfaces of 
different materials are superposed, meshes are 
20 repartitioned so that the shapes of the meshes of the 
surfaces are aligned. 

Fig. 35 exemplifies this mesh shape alignment. In 
his figure, a surface of a different material, which 
has a smaller area and is indicated by thick lines, is 
25 superposed on a surface, which is partitioned into three 
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meshes for respective pairs of opposite sides and 
indicated by thin solid lines. 

The mesh shape alignment is made by extending the 
edges of the surface indicated by the thick solid lines 
to the ends of the meshes of the large surface, and by 
repartitioning the meshes through which the extended 
lines pass in correspondence with the extended lines. 
As a result, the meshes before being repartitioned are 
deleted. For example, the upper right mesh is 
repartitioned into 6 meshes among the thin solid line 
meshes before being repartitioned in Fig. 35. Such a 
mesh shape alignment method is applicable, for example, 
to the case where a metal surface has a hole or a chink. 
In this case, the mesh having the hole or the chink can 
be handled as a mesh whose electric conductivity is 0, 
after being aligned. 

Once the process starts in Fig. 34, surfaces of 
different materials, which are superposed within a mesh 
area, are first extracted from mesh data and data of 
all surfaces in step S51. In step S52, it is determined 
whether or not the process is completed for all of such 
surfaces. If the process is not completed, an edge of 
the superposed surface is extended in step S53. Then, 
in step S54, meshes in the extended edge are 
repartitioned. In step S55, the repartitioned meshes 
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are added to the mesh data. In step S56, it is determined 
whether or not the process is completed for all of the 
edges. If the process is not completed, the operations 
in and after step S53 are repeated for the next edge. 
5 If it is determined that the process is completed 

for all of the edges of the extracted surface in step 
S56, the process goes back to step S51 where the next 
superposed surface is extracted. Then, the operations 
in and after step S52 are repeated. If the next 

10 superposed surface is not extracted in step S51, it is 
determined in step S52 that the process is completed 
for all of the surfaces. Here, the mesh shape alignment 
process is terminated. 

The above described generation of data input to 

15 an electromagnetic field intensity calculating device 
can be executed by a general computer, as a matter of 
course. Fig. 36 is a block diagram showing the 
configuration of such a computer system. In this figure, 
a computer 51 is configured by a main body 52, and a 

20 memory 53. As the memory 53, a storage device of various 
types, such as a random access memory (RAM) , a hard disk, 
a magnetic disk, etc. is available. A program recited 
in claims 7 and 8, or a program represented by the 
flowcharts shown in Figs. 14, 25, 32, and 34 is stored 

25 in such a memory, and the program is executed by the 
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main body 52, so that the generation of data input to 
an electromagnetic field intensity calculating device 
in the preferred embodiment can be implemented. Namely, 
it becomes possible to partition a metal surface of a 
5 metal cabinet of an electric circuit device into 
quadrilateral meshes which are suitable for a 
calculation of the moment method, and to output the mesh 
data after being partitioned to the electromagnetic 
field intensity calculating device. 

10 Such a program may be transmitted from a program 

provider side via a network 54, and executed by the 
computer 51. Additionally, such a program may be stored 
on a marketed and distributed portable storage medium 
55, so that the program stored on the portable storage 

15 medium 55 can be executed by being loaded into the 
computer 51 . As the portable storage medium 55, a storage 
medium of various types, such as a CD-ROM, a floppy disk, 
an optical disk, a magneto-optical disk, etc. is 
available. The above described program is stored onto 

20 such a storage medium, and executed by the computer 51, 
whereby the generation of data input to an 
electromagnetic field intensity calculating device, 
which is explained in the preferred embodiment, can be 
implemented. 

25 As described above in detail, according to the 
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present invention, a metal cabinet model that is 
suitable for a calculation of the moment method, namely, 
a model partitioned into quadrilateral meshes can be 
automatically and quickly generated from 
three-dimensional structure data of a CAD system. Such 
a model is output to an electromagnetic field intensity 
calculating device, thereby accurately calculating the 
intensity of an electromagnetic field emitted from the 
whole of an electric circuit device including the metal 
cabinet. This greatly contributes to improvements in 
the practicality of an electromagnetic field intensity 
calculating device. 



